clear all;

% 指定 Excel 文件路径和名称
filename1 = 'aerodynamic drag.xlsx';
% 指定要读取的工作表名称
sheet1 = 'throttle_speed';
% 指定要读取的数据范围，例如：'A1:B10'
dataRange1 = 'A2:B1869';
% 使用 xlsread 函数读取数据
sheet2 = 'throttle_lift';
% 指定要读取的数据范围，例如：'A1:B10'
dataRange2 = 'A2:B1869';
% 使用 xlsread 函数读取数据

sheet3 = 'speed_lift';
% 指定要读取的数据范围，例如：'A1:B10'
dataRange3 = 'A2:B1869';
% 使用 xlsread 函数读取数据

sheet4 = 'angle_lift';
% 指定要读取的数据范围，例如：'A1:B10'
dataRange4 = 'A2:B1869';
% 使用 xlsread 函数读取数据


[data, headers] = xlsread(filename1, sheet1, dataRange1);
% 将数据存储到 MATLAB 的变量中
throttle1 = data(:, 1); % 第一列数据
motor_speed = data(:, 2); % 第二列数据
% 油门转速


[data, headers] = xlsread(filename1, sheet2, dataRange2);
% 将数据存储到 MATLAB 的变量中
throttle2 = data(:, 1); % 第一列数据
lift = data(:, 2); % 第二列数据
% 油门升力

[data, headers] = xlsread(filename1, sheet3, dataRange3);
% 将数据存储到 MATLAB 的变量中
speedx = data(:, 1); % 第一列数据
lift_y = data(:, 2); % 第二列数据
% 转速升力



[data, headers] = xlsread(filename1, sheet4, dataRange4);
% 将数据存储到 MATLAB 的变量中
angle = data(:, 1); % 第一列数据
lift_ay = data(:, 2); % 第二列数据
% 转速升力

subplot(2,2,1);
plot(throttle1,motor_speed,'r-');
title("油门转速曲线");
xlabel("油门");
ylabel("转速");
legend('油门-转速曲线');

subplot(2,2,2);
plot(throttle2,lift,'r-');
title("油门升力曲线");
xlabel("油门");
ylabel("升力");
legend('油门-升力曲线');

[p,S] = polyfit(speedx,lift_y,2);
[fit_y,delta] = polyval(p,speedx,S);
subplot(2,2,3);
plot(speedx,lift_y,'ro',speedx,fit_y,'b-',speedx,delta,'g*');

title("转速升力曲线");
xlabel("转速");
ylabel("升力");
legend('原数据','拟合数据','拟合偏差');


subplot(2,2,4);
plot(angle,lift_ay,'r--');
title("(油门)舵角升力曲线");
xlabel("舵角");
ylabel("升力");
legend('(油门)舵角升力曲线');

grid on;


